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^vq ■ Abstract. We construct the quasi-classical approximation of the form factors in finite volume using the 

CO | separation of variables. The latter is closely related to the Baxter equation. 
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1 Introduction. 

00 

There is an important open problem of describing the matrix elements of local fields taken between two 
eigenstates of Hamiltonian (form factors ) for integrable field theory with periodical boundary conditions. 
There are at least two reasons why this problem is interesting. First, the knowledge of such matrix elements 
would allow us to study the correlators at finite temperature. Second, detailed understanding of the periodical 
problem would give controllable interpolation between a massive field theory and its ultra-violet, conformal, 
limit. 

In the present paper we shall study purely conformal case, more exactly, c < 1 models of Conformal Field 
Theory (CFT). Certainly, in CFT the correlators can be found explicitly, so one might find our study to be 
rather of academic interest. However this is not quite the case. The point is that we study the conformal 
field theory in its integrable formulation. The latter allows deformation to non-conformal case. 

The integrable structure of CFT was first discussed by Zamolodchikov [Q who constructed examples of 
higher local integrals. The existence of infinitely many local integrals was proven in ||. The spectrum of 
the local integrals is a subject of detailed study in the series of papers H 0]. In particular, the paper JD 
provides detailed investigation of famous Baxter equation. We are trying to combine the results of this 
paper with the method of separation of variables in integrable models developed by Sklyanin ||. In the 
latter method the solutions of Baxter equations play a role of wave functions for separated variables. We 
also use intensively the relation to the classical periodical solutions which are related to Riemann surfaces 
on which the separated variables represent divisors. 

The CFT compactified on the circle depends trivially on the length L of the circle. However, the case 
L — > oo is to be considered separately. Simple reasoning shows that the limit of the matrix elements in this 
limit must reproduce known form factors for Sine-Gordon (more exactly restricted Sine-Gordon) model in 
infinite volume Q. Here we find a relation to the paper || where the formulae for the form factor in the 
infinite volume have been explained in terms of separated variables. It should be said, however, that in the 
paper @] we failed to reproduce the quasi-classical limit of form-factors completely. More careful approach of 
the present paper will allow to find the missed pieces which are, in fact, due to the contribution of "vacuum" 
particles. 
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As it will be clear from the paper, though our results are quasi-classical, exact quantum formulae are to 
be found by similar means. How efficient these formulae would be is another question. 

Finally, I would like to say that the mathematical machinery used for study of periodical solution is the 
theory of Riemann surfaces. In L — > oo case these surfaces turn into degenerate surfaces which correspond 
to soliton solution. What should we achieve by studying the quantum periodical problem is a certain 
deformation of Riemann surfaces. For soliton case this deformation is explained in 0], [|l0); we hope to have 
something even more interesting in periodical case. 



2 The formulation of the problem. 

In this paper we shall consider the c < 1 CFT compactified on the circle of length L. We have the Virasoro 
algebra with generators C n satisfying the commutation relations: 

[C m ,C n ] = (m - n)C m + n + c 5 mt . 



12 
where 

c=13-6[S + r 



f K 

The Virasoro generators are obtained by usual construction from the Bose field <p(x) defined as follows 

<p{x) =iQ+ — — P + > — e — 

where the generators of the Heisenberg algebra satisfy the commutation relations: 

2,1 2 

The field ip is quasi-periodic: 

ip(x + L) = ip(x) + iLP 
This Heisenberg algebra allows the representsation with highest vector \p) : 

a n \p) = 0, n > 0; P \p) = p \p) 

later we shall use another notation for the zero-mode: 4> = —p. 
From this Bose field we construct the operator T(x) : 

T(x) = l(-. V >(xf:+(l-h)< f "(x) + ^ 

The representation of the Virasoro algebra in the space of representation of the Heisenberg algebra is defined 
because the Fourier components of T(x) satisfy the Virasoro commutation relations: 

\n= — 00 / 

obvioulsly, T(x) is periodical function of x. 

As it is shown in jl| the CFT allows integrable structure. It means that in proper completion of the 
universal enveloping algebra of Virasoro algebra one can find a commutative subalgbra of local integrals of 
motion. Namely, there exist an infinite sequence of local operators T2k{x) (for k = 1, 2, • • •) such that the 
operators 

L 

hk-i = 


commute with each other. The local operators T^x) are from the module of 1, i.e. they are costructed by 
taking derivatives of T{x), multipying and normal ordering them. In particular 

T 2 (x)=T(x), T 4 (x) =: T 2 (x) : 

c+2 
12 



J T 2k (x)da 



T 6 (x) =:T 3 (x):+^:(T'(x)) 2 : 
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The first local integral (Ji) coincides with C$. The spectrum of £q is highly degenerate, but other local 
integrals of motion reduce drustically this degeneracy. There are two important problems. First, one has to 
describe this spectrum. This problem is very complicated, but as it is conjectured in the spectrum can 
be defined from solutions of Baxter equations (we shall discuss this later). 

The second interesting problem consisits in calculating the matrix elements of local operators between 
the eigen-states. Suppose that we have two eigen-states \^) and \^f') such that 

/ afc -i|*) - iafc-il*), / afc -i|*') - iifc-il*') (1) 
where i, i' are eigen-values. We are interested in the matrix elements of the kind 

(9\0(x)\V) 

where 0(x) is certain local operator. For simplicity we can take 0(x) from the module of 1, in that case, 
obviously, the states must have equal zero-modes, otherwise the matrix element vanishes. In this paper we 
shall consider the matrix elements between two states with <f> = n. Obviously, it is enough to consider the 
matrix elements of O(0) because |$) and are eigen-states of Cq, so, the x dependence is trivial. 

What can be said in general about the matrix elements ([!])? There are two situations when they are 
known. 

1. Consider the "small" states, i.e. the ones created from Virasoro vacuum by applying few rising operators. 
On these states the local integrals can be diagonalized explicitely, and the matrix elements can be found by 
brut force. 

2. More interesting case is the case of "big" states. There are two equivalent definitions of these states. It 
is a proper place to explain why we keep the lenght L in all formulae. In conformal theory we can always 
rescale the circle to one of length 2ir passing to the variable y=^j-. The local integrals must be rescaled as 
follows 

L ^ 2k -' 



Ja *- 1 - hh - x 

where I are the integrals on the 27r-circle. Consider the the states on i-circle for which the eigen-values of 
the local integrals remail finite in the limit L — > oo or, equivalently the states on 27r-circle for which the 
eigen-values of I 2 k-i are of order L 2fc_1 for certain big parameter L. This is the definition of "big" states. 
We prefer the first interpretation that is why we keep L in our formulae. It is rather clear that in L — > oo 
limit the matrix elemnts between "big" states must reproduce the form factors in restricted SG theory. 
Indeed the large L limit corresponds to conformal thory in infinite volume. This theory can be described by 
massless S-matrices |5j and form factors which don't differ from those of massive theory. The states with 
<f> = n correspond to n-soliton states in infinite volume limit. 



3 The classical periodical problem for KdV. 

Let us present several facts concerning the classical KdV hierarchy with periodical boundary conditions 
following mostly the book We have the field u(x) satisfying the periodicity conditions u(x + L) = u(x), 
the second (Magri) Poisson structure is defined by 

{u{x\), u(x 2 )} = (u(xi) + u(x 2 ))S'(xi - x 2 ) + e5"'(xi - x 2 ) (2) 

We shall consider two "real forms": the case e = + corresponding to usual real solutions of KdV equation 
(this case is denoted by rKdV), and the case e = — which corresponds to certain class of complex solutions 
(this case will be denoted by cKdV). The latter case is related to the c < 1 models of CFT discussed in 
the previous section because in the quasi-classical limit the Poisson brackets of u[x) = TiT(x) coincide with 
. Actually the two cases must be understood as analytical continuations of each other, see [|| for more 
expanations. 

The KdV hierarchy possesses infinitely many integrals of motion I 2 k-i (k > 1) the first of them (the 
momentum) being 

L 

h =Ju(x)dx 
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Making the rescaling 

y = "(y) = l2u (J^y 

we map the periodical problem for u with arbitrary L to the one for u with the period equal 2-7T, the integrals 
of motion scale as in the quantum case. 

The exact solution of KdV equation is due to existence of Lax representation. The auxiliary linear 
problem is 

^-u(z))<Mz,A) = A 2 ^(x,A) (3) 

It is convenient to rewrite this equation as matrix first order equation. To this end introduce the field 
related to u(x) by Miura transformation: u — (<p') 2 + f" ■ The field ip is real for rKdV and pure imaginary 
for cKdV. The equation (||) is equivalent to the following one: 

|*(i,A)=£(i,Mi ) A) 
ax 



where 



V e^r 1 (ip(x, X)' - ip(x, X)cp(x)') 



and 

£(x, A) = 



P -¥>(a0 

2 ' 



\2 

, c ' 2 

The fundamental role is played by the monodromy matrix: 



(x +L \ 
J C(x,X)dx\ 



The trace of T XQ (X) which does not depend on xq is denoted by T(A). Let us recall the main properties of 
T(A). The function T(A) is an entire function of A 2 with infinitely many zeros accumulating to oo along the 
negative real axis. The simplest example is given by u = for which 



T(A) = 2cos(iV^A2) 
The graph of this function looks as follows: 




Figure 1 

Further we shall give other examples. The monodromy matrix is a 2 x 2 matrix with unit determinant, 

so, 

T(A) = A(A)+A(A)- 1 

where A, A -1 are the eigenvalues of T. It is important to notice that A(A) is not an entire function of A, it 
has quadratic branch points at zeros of odd order of the discriminant A(A 2 ) = T(A) 2 — 4. The discriminant 
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is an entire function of A 2 with asymptotic following from (^J), so, it can be described by converging infinite 
product over its zeros 

where C = dA/d\ 2 (0), kj = 1,2. The asymptotical behaviour of A(A) is governed by the local integrals of 
motion: 

logA(A)~ \L + Y,^ 2k+1 hk-i 

k>l 

A -► oo, RcA > (4) 

The function A(A) is a singe- valued function on the hyper-elliptic Riemann surface £ (generally of infinite 
genus) given by the equation 

= a 2 n (a 2 - «*> 

fci=i 

The tractable mathematically and (fortunately) at the same time most interesting physically case is when 
£ has finite genus, i.e. when A(A 2 ) has only finitely many zeros of first order. The typical situation of this 
kind is given by the periodical analogues of n-soliton solutions. It that case we have simple zero at A = 0, 
2n real positive single zeros of A (we denote them by A 2 , • • • , \\ n and infinitely many negative double zeros 
(we denote them by — /J?_i, — M-2 ' ' ')• The graph of the function T(A) looks as follows: 



Figure 2 

The segments of the real axis of A 2 where |T(A)| > 2 are called forbidden zones of the periodic potential 
u (there are no bounded wave- functions for these values of energy). 
The hyper-elliptic Riemann surface E of genus n is described by 

2n 

/i 2 = A 2 F(A), P(X) = l[(X 2 A 2 ) 

4=1 

Conventionally it is realized as two-sheet covering of the A 2 plane with cuts and canonical basis of homology 
chosen as follows 




■5 
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Figure 3 

We shall be mostly using another model of this surface considering it as A-plane with cuts: 

in 



in 



in 



-X, 4 —X, 3 —A 2 — A i 



a 

i l 




Figure 4 

The upper (lower) bank of the cut [Aoj—i, Ajj] is identified with the upper (lower bank of the cut 
[— A2J-1, — A2j]. Obviously the upper (lower) half planes correspond to first (second) sheet of the surface in 
usual realization. 

The function A(A) continues to the lower half-plane as 

A(-A) = A(A)- 1 

Let us consider in some details properties of this function. 

The function log A(A) is called quasi-momentum because the equation (^) has the Floquet solution ij>(x, A) 
which satisfies quasi-periodicity condition: 

ip(x + L,\) = A(X)ip{x,X) 

another name for this solution is Baker- Akhiezer (BA) function. The BA function is single-valued function 
on S, ip(x, A) and ip(x, — A) give two linearly independent for generic A quasi-periodical solutions of ([|), the 
second one corresponds to the quasi-momentum — log A(A). 

The function A(A) is single- valued on £ with essential singularity at infinity. Hence dlogA(A) is abelian 
differential on £ with second order pole at infinity: 

dlogA(A) = (L + 0(A" 1 ))dA 
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(the local parameter at oo is A 1 ). Moreover it must be normalized: 

dlogA^O Vj (5) 



because as it is seen on fig. 2 A(A) is real positive or negative function in the interval corresponding to a-cycle. 
Then in order that A(A) is single-valued we need that 

d\ogk = 2mk j (6) 



for some integer kj. This requirement can not be satisfied for arbitrary S, it is actually a restriction on 
the moduli of the surface. For the periodical analogues of n-soliton solutions we have kj = n — j + 1 which 
corresponds to the fact that exactly n — j + 1 simple zeros of T(A) are surrounded by bj. This is exactly the 
situation presented on the fig. 2, in more general case between two forbidden zones T(A) can make several 
oscillations from -2 to 2. Provided and (Q) are satisfied logA(A) is a function defined on £ with cuts 
along the a-cycles whose jumps on aj equal 2nikj. 

The dynamics of finite-zone solution is conveniently described by motion of zeros of BA function. 
There are exactly n of them (70, • • • , 7„_i). Let us present for completeness the equation describing the 
x-dependence of 7,*: 



The dynamics with respect to higher times is described by similar equations. The dynamics is linearized by 
Abel transformation of the divisor of zeros of BA function onto Jacobi variety of E. Recall that we consider 
two different real forms of KdV (rKdV and cKdV) . The points of divisor corresponding to this two real forms 
move along topologically equivalent, but geometrically different trajectories. In rKdV case jj moves along 
the cycle aj as it is drawn on the fig. 4. In cKdV 7; runs along a trajectory close to the cycle Cj on the fig. 4. 
Clearly, the half-basis of c-cycles is equivalent to the half-basis of a-cycles. 



4 Hamiltonian structure of finite-zone solutions. 

Let us discuss the most subtle issue in the theory of finite zone integration, namely, the Hamiltonian de- 
scription of the solutions. The surface E is parametrized by In parameters (Ai, • • • , A2 n ). These parameters 
are not all independent, they are subject to n restrictions (j^. So, we are left with n-dimensional sub- 
manifold {Ai) in the moduli space of hyper-elliptic surfaces. It is convenient to parametrize M. by ti, • • • , r n 
such that t 2 are positive zeros ot T(A) in A 2 -plane. Earlier we have introduced the variables 70, • • • , 7n-i- 
The phase space of finite-zone solution is the 2n-dimcnsional manifold locally described by the coordinates 
{ti, • • • , t„, 70, • • • , 7n-i} is embedded into the infinite-dimensional phase space of KdV. Restricting the sym- 
plectic form which corresponds to Poisson structure (|^) to this finite-dimensional manifold we obtain the 
symplectic form ui = da with 1-form ui given by 

n — 1 , 

a = 5>gA( 7j )^L 

Our main concern is quantization, so, we have to ask ourselves the question whether the quantization of 
this finite-dimensional mechanics is relevant to the real quantization of KdV. Logically, the answer to this 
question is negative because restricting ourselves to the finite-dimensional sub-manifold we ignore a good 
deal of quantum fluctuations allowed in the infinite-dimensional phase space. However, as is shown in | flo| 
for the case of solitons (L — > 00 limit) the quantization of the finite-dimensional system gives an important 
piece of exact quantum answer for the matrix elements (form factors 0). To understand why it works and 
how to generalize the results of [jioj] to periodical case we have to consider the hamiltonian structure in some 
more details. This consideration will allow also to reproduce quasi-classically an important part of solitons 
form factors which we could not do in fllOfl. 
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Let us analyze more general situation of which KdV provides a particular case. Take a class of clas- 
sical integrable models with trigonometric R-matrix. For any such system (continuous or lattice one) the 
monodromy matrix 

can be introduced which satisfies the famous Sklyanin's relations 

{T(A) f T( M )} = [r(A, a*), T(A) ® T(a*)] (8) 

with trigonometric R-matrix: 





For a lattice regularization of KdV |T2| the monodromy matrix is a polynomial in A 2 of degree N. The deter- 
minant D(\) of T(A) is in the center of the Poisson algebra. The trace T(A) of T(A) is a generating function 
of jV independent integrals of motion. To describe the phase space one has to introduce N coordinate-like 
variables. Following general approach developed in jl^, |(| we take as such zeros of B(X): 

N 

B{\) = B(0) n 

3=1 

Using ( |TT| ) one finds: 

{A( 7i ), A( 7 ,)} = { 7l , 7j } = {A( 7i ), 73} = 5 lj7l A( 7l ) 
where A(A) is eigenvalue of T(A): 

A(A) = IW±^M 

the discriminant A(A) = T(\) 2 — 4D(A). Thus the variables log 7 j and logA( 7 j) are canonically conjugated 
and the symplectic form can be written as lo = da with 



N A 



a = ^logA( 7i )-^ (9) 

3=1 7j 

Generally, the dynamics is linearized on the Jacobi variety of the hyper-elliptic Riemann surface E described 
by /j, 2 — A(A). The genus of E equals N, the points of the divisor 71, • • • , 7 at move along certain closed 
curves topologically equivalent to a half-basis of a-cycles on E. 

Later we shall need also the Liouville measure corresponding to this symplectic form. Taking zeros of 
T(A) (tjl, • • • , tjv) and 71, ■ • • , Jn as coordinates on the phase space one finds 

aN " = U a ( \ \-K s li^-^Ii^-^— ^■■■^^^drl^...^dr N (10) 
fJl Mw - a {wt^j f<j 71 1N 

What happens to all that in the continuous limit which corresponds to finite-zone solution of KdV? The 
degree of T(X) goes to infinity, but the surface E turn into a surface of infinite genus of rather special type. 
Namely, only finitely many of zeros of A(A 2 ) (0, A 2 , • • • , A 2 n ) remain simple ones while infinitely many zeros 
— /i 2 _2i ' ' ') become double zeros. The polynomial T(A) looks on the lattice as 
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Figure 5 



In the continuous limit it turns into the one presented on the fig. 2. In the lattice regularization of 
rKdV case the points 7,- are moving inside the zones where |T(A)| > 2. In N — > oo limit only n of them 
(70) • • • , 7n-i) move inside the |T(A)| > 2 zones on A 2 > half-axis on fig. 2, while infinitely many (we shall 
denote them by 7-1,7-2, ■ • •) happen to be confined at the points — /r?^. This is how the restriction of degrees 
of freedom takes place in the classical case. It must be emphasized that this is a dynamical procedure which 
can not be carried out on quantum level. The points 7-1, 7-2, • ■ • can not be kept at fixed positions, but 
must be distributed with certain wave function localized around their classical values. 

To finish this section let us present some explicit formulae concerning zeros of B(X) in the continuous 
case. Consider the equation (||) with finite-zone potential u{x). Let us fix the normalization of the BA 
ip(x, A). If we require that 

ip(X) ~ exp(Aa;), A — * 00 

the BA function must have n simple poles in the finite part of the plane. It is convenient to put this poles 
to n branch points, say, Ai, • • • , A n . We can introduce another function ip^x, A) which satisfies 

^(x + L,X) = A- 1 (A)V t (a;,A), V f (A) ~ exp(-A^), A -> 00 

and has poles at the complimentary set of the branch points. These two solutions satisfy the relations 

nV-7?) 

W(^(X),^(X)) = 2X, V(A)V f (A) = 3 

where W(tp, ift) is the Wronskian determinant. From these two solutions to (||) one easily reconstructs the 
monodromy matrix finding in particular that 

B(X) = (A(A) - A-^A)) ^(A)^(A) - B(0) JJ l_^ JJ i + 4_ 

j=0 \ ~j ) j=l \ V-j j 

Thus we find that 6(A) has indeed zeros not only at the moving points of the divisor, but also at the confined 
points fi-j. 



5 Separation of variables. Baxter equation. 

Our nearest goal is to write down the quasi-classical expression for the wave function corresponding to the 
periodical analogue of 71-soliton solution. As we have seen there are two types of coordinate- like variables: 
7i, ' ' ' , 7« which moves classically along the a-cycles and 7-1,7-2, • ■ ■ which are classically confined. The 
first type of variables does not pose a problem for writing the quasi-classical wave-function. The contribution 
to the wave-function from the second type of variables is similar to that of a number of harmonic oscillators 
in the ground state. To understand this contribution we will need some pieces of exact quantum information. 

Consider an operator-valued monodromy matrix T(A) with the same notations for the matrix elements 
as before which satisfies the quantum analogue of the Poisson brackets (||): 

R(X, /i)(T(A) <g> /)(/ ® T{jj)) = (I <g> T[f Jl ))(T(X) ® I)R(X, p) (11) 
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where 

/X 2 q-/i 2 q- 1 0, 0, 

0, A 2 -^ 2 , ^(q-q- 1 ), 

H ^^>- 0, \ 2 {q-q- r ), A 2 -^, 

V 0, 0, 0, X 2 q-fl 2 q~ 1 

where 

q = e** R , 

h is the Plank constant. In conventional notations coming from SG theory 

» = A 

The quantum determinant of T(A) is in the center of the algebra, in quantum KdV case it is fixed to be 1: 

A(X)V(qX) - B{X)C(qX) = 1 

where q — exp(i£). The monodromy matrix T(A) is an entire function of A 2 . The trace of the monodromy 
matrix T(A) is a generating function of integrals of motion. In quantum KdV case one can construct the 
monodromy matrix directly in continous theory j^] by proper normal ordering of the classical formula. 

Let us review the method of separation of variables which was developed by Sklyanin raj who combined 
the general approach of Inverse Scattering Method H] with the ideas of and [H). In application to 
our particular case the results of Sklyanin can be formulated as follows. Consider the element 23(A) of the 
monodromy matrix. It defines a commutative family of operators due to the fact that 

[23(A),23(/i)] = 

Moreover 23(A) is an entire function of A 2 which grows not faster than exp(A p L) (with certai p, see below) 
at the infinity, that is why it can be presented as an infinite product over its zeros: 

B(A) = 23(0)fj (l- (^) J (12) 

The idea of Sklyanin is to consider the system in "7-representation. 

The functions -4(A) and T>(X) are entire functions of A 2 and as such they can be expanded into series of 
infinite radius of convergence: 



A(X) = J2 >? n An, V{X) = A 2 "2?„ 



n=0 n=0 

Following H one introduces the operators 

A j = A(j j ), Aj = V(7j) 

where the operator jj is substituted into A and T> from the left (exactly replacing A in ([L2|)). The operators 
7, A, A posses the following important properties: 

[7j,7fc] = 0, [Aj, A fe ] = 0, A i7fc = q^-fkAj, AjAj = 1 

first three equalities follow directly from the commutation relations while the last one is the consequence of 
the commutation relations together with the quantum determinant. These commutation relations show that 

( ■ * 9 \ A ( ■ * 9 

Aj = ej exp tirn~fj— — , A., = e,- exp —nrnr/j— — 
V 3 &YjJ V dlj 

where tj = ±1, the necessity of introducing tj comes from consideration of the fig. 2: the variables 7 move 
(or rest) classically inside zones where |T(A)| > 2, so Cj = ±1 depending on whether T(A) > 2 or T(A) < —2. 

Our goal is to diagonalize the Hamiltonians (the operator T(A)) in the 7-representation. This can be 
done if we accept the following: 
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Conjecture. 

The function B(X)^ 1 T(X) is a meromorphic function of X 2 with infinitely many poles which can be expanded 
in convergent series: 

n— 1 

B(X)- 1 T(X)= ]T -— ^(A.+A,) 

j=-oo A '3 



Assuming that the wave-function is presented as the infinite product: 

n-1 
j=-oo 

where every Qj satisfies the equation: 

e j t(X)Q j (X) = Q j (qX)+Q j (q- 1 X) 

The sign ej can be taken away multiplying Qj by appropriate power of A, and basically we have to study 
the famous Baxter's equation: 

t(X)Q(X) = Q(qX) + Qig-'X) (14) 
where t(X) is the eigenvalue corresponding to \I/: 

T(A)¥ =t(A)« 

All this procedure is called the separation of variable because it allows to transform the infinite-dimensional 
spectral problem to one-dimensional ones. The separation of variables in quantum case gives an exact 
quantum analogue of the classical separation of variables which is obvious in the formula (|^) . 

The equation fll4|) can be considered as a finite-difference analogue of Schrodenger equation, this analogy, 
however, is not quite straightforward. It the case of usual Schrodenger equation the wave-function belongs 
to certain functional space like L%. In our case the wave-function Q(A) is characterizes by its analytical 
properties. 

As a second order difference equation ( |l4| ) must have two solutions up to multiplication by quasi-constant 
(a function of An). Following Q we require that one of them (denoted later by 2(A)) is an entire function of 
A 2 . Recall that t(X) is also supposed to be an entire function of A 2 . The next important piece of information 
is the asymptotic behaviour of these functions. The equation ( |l4| ) can be rewritten as 

t(X)=A q (qix) + A-\q-ix) 

where 

A 9 (A) = (15) 

is the quantum analogue of the eigenvalue of monodromy matrix. In quantum KdV theory A g (A) allows the 
following asymptotical expansion around infinity: 

log A,(A) ~ X P L + i<f>ir + X- {ik+1)p hk-i 

k>l 

X — > oo, nh < argA 2 < n(l - K) (16) 

where 

1 

P 



i-h' 

^2fc-i are quantum local integrals of motion (see for exact normalization of them). The fact that T(A) 
is an entire function of A 2 that has the asymptotical expansion in terms of X p is well known [|l4). The 
explanation of this fact is due to the renormalization of mass. More modern and clear explanation in terms 
of CFT is given in B. Notice that our definition of A coincides with the one used in B. The phase <j) 
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equals to the topological charge, for the periodic analogue of n-soliton state <fi — n. The leading terms of the 
asymptotics of Q(A) are supposed to be: 

logQ(A) — j\ p L—\cf) logA + (3(A~ p ), A 2 ^oo, < argA 2 < 2tt (17) 

2i sin | n 

In the paper M one finds detailed description of complete asymptotical series for Q. 

Furthermore it is required that the function i(A) has only finitely many zeros away from the axis A 2 <0 
(n real positive zeros for the periodic analogue of n-soliton state) while the function Q(A) has only finitely 
many zeros away from the axis A 2 > (for the periodic analogue of rt-soliton state they are absent). The 
basic conjecture accepted in j| is that counting the solutions of ([li]) one counts the vectors in the space of 
states of CFT. This conjecture is verified in many cases, so we have no doubt that it is true. 

The entire function Q(A) can be presented as infinite product with respect to its zeros: 

Q(A) = Q(0)n^l-(A^ (18) 

The zeros Oj are subject to the Bethe Ansatz equations 

Qiqa ^ = 1 (19) 



In the next section we shall consider the solutions to these equations in the quasi-classical limit. 

As it has been said there must be another solution (Q(A)) to the equations (pT[). For any two solutions 
to (|l4|) the "quantum Wronskian" 

W(Q, Q)(A) = Q(X)Q(Xq) - Q{Xq)Q(X) 

is a quasi-constant: 

W(Q,Q)(qX)-W(Q,Q)(X)=0 

So, to find the second solution we have to solve a first-order difference equation. Namely, let us put 
W(Q, Q)(A) = 1 then Q(A) = Q(A)F(A) where F satisfies the equation 

FW - F{X> = mkm (20) 

The solution to this equation can be always found but generally it is not a single- valued function of A 2 . 
Consider for simplicity the reflectionless case Ti = i, with integer v. It this case Q is a single- valued function 
of A 2 described as follows. Take some polynomial P{a) of degree n then 



F(A) - P{X2V) hi J Q^QWW^^-A 2 ) + hi J PW) g(Ai m) 

\ Ci c 2 / 

where C\ encloses zeros of Q([i) and C2 encloses zeros of P(/j, 2l/ ), G(A, /z) is the following function: 



(21) 



v 3 



G(A,M) = EE 



^Z_. (aVj -_ m2)q(m -; )q(m -; +1 ) 

The polynomial P is introduced for convergence, we shall not go into more details here. The definition of F 
depends on P, but one easily shows that solutions with different P differ by a quasi-constant. 

6 Quasi-classical wave functions. 

Our nearest goal is to understand the quasi-classical behaviour of Q(A). Consider the equations: 

_! (22) 
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In the case which we shall consider in this paper all the solutions to these equations are such that A 2 is real. 
A part of solutions coincides with zeros of Q(A) but, obviously there are other solutions which provide zeros 
of t(X). Different solutions are counted as follows: 

log ( - + 1)- 



where — oo < k < oo. We have to share these solutions between t(X) and Q(A). For the quantization of 
periodical analogues of soliton solutions we do it as follows. For — oo < k < Xk are zeros of t(A), for 
< k < oo they are zeros of Q(A) except for finitely many N\, ■ ■ ■ , N n > which correspond to positive 
zeros of t(X). We have taken the branch of logarithm such that the border between zeros of t and zeros of 
Q lies at k = 0. 

In quasi-classical limit K — > 0, Nj = Olfi -1 ) and the zeros of <2(A) condense in the A 2 -plane forming the 
cuts |l6| . From comparing with classical picture it is clear that these cuts must coincide with the intervals: 
I\ = [0, A 2 ], I2 — [A§, A§], • • • , I n +i — [A 2 „, 00] where A 2 are the branch points defining the classical solution. 
There is one zero of t(X) in every interval between Ij and Ij+i- The quasi-momentum logA(A) can be 
considered as single-valued function on the plane with these cuts. Comparing classical formula (^|) and 
quantum formula (Q3|) one easily finds that when q — > 1 



A 

logQ(A) = ±- ( \ogA(a)—+ O(h ) 
inn J a 



This formula must be understood as very approximate one, logA(A) is not a single- valued function, and its 
branch has to be taken differently for different jj. This "tree approximation" of Q(A) can be rewritten in 
the form: 

x 7 



where inside every Ij the density 



P (a) = -Re(logA(a)), 
a 



notice that p(a) > inside Ij, and it vanishes as yj\a — Xj \ at the ends of intervals. The logarithms in ( J25 
have cuts in a 2 plane from A 2 to 00. Obviously the formula (|2^) originates from the classical limit of the 
infinite product (|Tg|), and p(a) describes the density of distribution of zeros. 

We need the quasi-classical correction to Q(A). To find it we have to investigate carefully the classical 
limit of Bethe Anzatz equations ([jl]). It is done in the Appendix using certain version of Destry-de Vega 
equations [0, here we present the final result of these calculations: 

Q(A) - \-ip~i(\) r s (A)exp J logA^^j (24) 

where Fs(A) is an analogue of gamma- function related to the surface S which is defined as follows. 
Let us introduce the normalized third kind differential /o^(A) = G^(X)dX 

dX f 

Then 

dlogr s (A) = - Jim ( -L / p,(X)d\og(A(p) 2 - 1) - Alog^(A) I (25) 

\ Cjv / 

where the integral is taken over /i, and the contour Cn encloses the points — ijij with 1 < j < N, 

p(X) = lim M ^ 00 /ip M (A) 
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The function Ts(A) is not single- valued on S, but it is single- valued on the plane with cuts (fig. 3). 

The analogy with usual T- function is obvious. If we take, in particular, the vacuum state the correspond- 
ing Riemann surface S uoc is Riemann sphere and 



LX 



The function satisfies the important functional equation: 

r s (A)r E (-A)- A(A) _ A _ 1(A) (26) 

Before going further into study of the wave functions let us consider the quasi-classical analogue of the 
quasi-momentum. By definition 

VA) = §^l v (27) 
Q(Xq 2) 

its asymptotics are (plf): 

log A,(A) ~ X P L + ^ (2k+1)p hk-i 
fc>i 

recall that p — j^r . Obviously one has quasi-classically: 

logA 9 (A) ~ R ^ logA(A) + fi ( log A A^ log A(A) + ]T X- (2k+1 ^5I 2k ^ ] (28) 

fc>i 

where Sl^k—i ar e quasi-classical corrections to the integrals of motion. The non-trivial term of this formula 
is that containing log A. Let us show that our previous formulae agree with this asymptotic behaviour. From 
one finds: 

log A,(A) = log A(A) + ihX^- log(A-iP(A)-ir s (A)) + 0(h 2 ) (29) 

To compare the formulae (|29| ) and (|28[ ) we need to know the asymptotics of Tx;. It is quite clear from the 
functional equation (^6|) that asymptotically 

1 oc 

log(A-4p(A)-ir s (A)) = -logA logA(A) + Vc fc A- 2fe+1 

fc=l 

where Ck are certain coefficients. Now we see that the formula (^) has indeed the same structure as 
Actually, the coefficient can be evaluated explicitly providing quasi-classical corrections to the integrals 
of motion. 

Finally, we have to require that Q(A) is single- valued on the plane of A 2 with cuts (fig. 3). This requirement 
leads to the Bohr- Sommerf eld quantization conditions: 

J dlogQ(A) = 2mNj (30) 

aj 

obviously Nj is the number of zeros of exact quantum Q in corresponding interval, in the quasi-classical 
region Nj are of order fi~ . These quantization conditions are implicitly the quantization conditions on the 
moduli t% , ■ ■ ■ , T n . 

The following important circumstance must be emphasized. The function Q constructed in pi is defined 
in A 2 plane, so, it has as its quasi-classical limit the function constructed above defined in A 2 plane or, 
equivalently, in the upper A plane. However in our construction this function allows analytical continuation 
to the lower half-plane (second sheet). This analytical continuation is related to the second solution of the 
Baxter equation Q discussed in the previous section. 
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Let us construct the quasi-classical wave function. We had the formula 

n-l 

*= n 

3=— oo 

As it has been said the functions Qj{X) differ from Q(A) by certain degree of A. Let us explain the origin of 
this difference. First, notice that the expression for Q (E4h is not uniquely defined on £ because d\og A has 
non-zero 6-periods (all of these periods are equal to 2ni, so adding them we would multiply <2(A) by At). 
Let us fix the branch of log A(A) as follows: make cuts over the cycles Cj (fig. 4) and require log A(0) = at 
infinity. Now if we understand the integral in (|24| ) as 

A A 

o 

The functions Qj(X) are defined as 

Q j (\) = \tQ(\) (31) 
This prescription is chosen for the following reasons. 

1. For j < it ensures the existence of saddle point of Qj at ifj,j, the point where jj is situated in classics. 

2. For j > it makes the action to satisfy proper reality condition along a classical trajectory. The problems 
of reality are discussed in details in || . 

The prescription for quantization of cKdV explained in [0 consists in taking the analytical continuation 
of rKdV. In particular it corresponds to the following rule of hermitian conjugation: 

Hrt( 7 )=-^7) (32) 

which is the same as 

where 7 C = 7 for classically confined coordinates and 7 C = —7 for classically moving particles. 



7 Quasi-classical matrix elements. 

To construct the matrix elements we need to know also the measure of integration in the space of functions 
of 7. Comparing the expression for the Liouville measure (0) with the formula @ we assume that the 
measure of integration is given quasi-classically by 

n-i 

w( 7 )= n ^-^) 

— oo<i<j 

Thus the matrix elements of operators in 7-rcpresentation are given by: 

n-l 00 

mom =^I| / d ^ W/(7)*(7 C )0(7)*'(7) (33) 

where JV, N' are norms of the wave functions, 0(7) is the operator O in 7-representation. Let us first 
consider the norms. 

Take the wave function for the periodical n-soliton solution ^(7) and consider 

n-l 



N2 = II / d ^ ^(7)*(7 C )*(7) 



Consider first the integral over 7-j. It is of form 

7 

da 



J F(j)exp I j (logA(o-) -j-ni)-^ 
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where F(X) is finite when h — > 0. So, it is sitting on the stationary point ifj,j (recall that log A(/ij) = jm). 

Now consider the integral with respect to jj (j > 0). By definition of branch of log A(A) one sees that 
on the real axis 

logA(A) + logA(-A) = ™fc for |A fc | < |A| < |A fe+1 | 



where we put Ao = 0, A2 n +i = oo. Together with (26) it gives 

^ g (-^ ¥ - A (7 )A-H 7 ) n/i 

for |A2j| < |A| < |A2j+i| and exponentially (with ^ in exponent) smaller everywhere else. So, the integral 
with respect to 7j is sitting on these two segment of real axis or, putting it differently, on the cycle aj. 

After these expansion let us present the final result of calculations. We denote by A q the quasi-classical 
approximation given by (29). 



2n £ n — 1 oo y 2 \ 

a' 2 = *(a 9 ) n xi 1 1 10 ■ ■ ■ 1 dln ^ n A(7j) j^ (7j) n (i - £) n w - 

k=l ao a„_i 3=0 k=l x " ' i<] 



2)1 



$(A ? ) n A (34) 



fc=i 



where 



and 



A = 



fa- J ^nJnW*-* 



*(A,) = exp ( ^ flogCA^) - 1) (£logA») lo gA1 + ^log (jfclogA^) 



d/i— 



-^//log(A2( Ml ) - l)log(A2(^) - iJ^WMMa) , (35) 

cc 1 2 / 

the contour C encloses the points ip-u- The formula ([m|) must be interpreted as follows: A gives the 
volume of the coordinate space of the finite-zone solution while the rest describes the measure in orthogonal, 
momentum, direction. Actually the integrals in (|35| ) are divergent, one has to divide by the norm of vacuum 
to make them finite. This divergence does not affect our further calculations, so, we shall ignore it. 

Let us consider now the matrix elements. We shall take the simplest operator T(0) which is classically 
the same as it(0). On a classical finite-zone solution u = Jj — \ S Aj • The prescription for the symbol of 
this operator in quasi-classical approximation is: 

where A, A' are branch points corresponding to classical solutions. We take this symmetric prescription 
because it looks the most natural and gives correct answer in L — > oo limit || . The calculation of the matrix 
element is similar to the calculation of norm. The only important point to realize is that the stationary 
points move to those solving the equation 

A(M)A'(/i) = 1 

The final result is 



OC OC / 

X / <*y -" / dTn-inQ^OQ'CT.OexpU/^MA^A'^-l)^ 

— oo — oo j \ C j 

x nil n (if - 1]) (Et| - 1 e( a2 + (a;-) 2 )) (36) 

where the contour C encloses the zeros of A(/x)A'(/x) — 1 lying on positive imaginary half-axis. Notice that 
the exponential growth at infinity is cancelled in the combination Q(— Jj)Q'(jj)- Before discussing further 
this formula let us show that it gives correct result in the limit L — > oo. 
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In the classical case when L — > oo the Riemann surface degenerates: with exponential in L precision 
Aaj-i — * Tj <— The quasi-momentum becomes an elementary function: 



A(A) = e iA 
In the quantum case one has @ : 



f[(|^)(l + 0(e-^)) (37) 



A 9 (A)=e- P n Mrp£ (l + 0(e-^)) (38) 



Let us first check that our quasi-classical formula for Q agrees with this result. Consider the formula (^5|). 
On the degenerate surface the differential p turns into 

^(a) = (tx^)+^- 1 ))^ 

One finds that when L — > oo 

dlogr s (A) = -lim^oc p M (A)dlog(A( A1 ) 2 - 1) - -Llog^p(A)^ -> 

lim M ^co i ^ / ^logA^x) - LlogM^ dA (39) 

where we have integrated by parts and took into account that A(A) is exponentially big (small) in the right 
(left) half plane. Calculating the integral in ([39j) one gets: 

d (A-^P(A)-i lo g r s (A)) -> 1 flog A rflogA(A) +^ 2 J| ^ d\j 

Now we substitute this expression into the definition (|27|). Obviously, 



Xh ^log A A bg A(A) + £ 2 r l l0g A ^ - log A, (A) - log A(A) + 0{% 2 ) 

which proves the consistency of our quasi-classical formulae with exact quantum formula (p8|). 

The calculation of the matrix element in L — > oo limit is straightforward, but bulky. The main simplifi- 
cation is due to the fact that the integrals containing log(A 2 (^) — 1) or log(A(/i)A'(/x) — 1) can be evaluated 
as it has been done in the calculation of dlogrs(A). The only non-trivial part of the calculation is that 
concerning A. Recall that the differential <ilogA is a normalized second kind differential on the surface. 
Requiring that A is given by fllT^ ) in L — > oo limit we actually fix completely the rule of degeneration of 
the surface in this limit. In particular, one easily finds the limiting values of the normalized holomorphic 
differentials, and shows that 

A _^ \ T n 

i<j 

The appearance of L n is not surprising because in L — > oo limit the eigenstates are normilized with 5- 
functions. 

To make easier the comparison of the final result with known exact formulae [Q it is convenient to 
introduce usual rapidity notations: 

0j = ilogTj, ay = ~ log 7,7 
Then the final result of the calculation is: 

<*ir(o)i*') = L n p p> n c(a - Pi) n m - %) n cc* - # - m) 

i<j i<j i,j 

71 — 1 TL TL 

xf ai ---fda n n n ^ -/%-?) n 

i=0 j=0 j=0 

x n ef^- n+1 ^ n sinh(a, - a,) fe e<* - \ £ - § £ e#) (40) 
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where 



OO OO — 7TZ 




— oo —oG — iri 



the functions £ and ip are given (up to some constants) by 

which is exactly the quasi-classical limit of corresponding functions used in [Q . Finally, 

p = cx P (- £(2j - n - m) [] y^H^-^) ~ & + y) 

and similarly for P'. Notice that \P\ = 1. 

One can show that the formula (|l|) basically coincides with the quasi-classical limit of modest modi- 
fication (similar to one done in ||) of usual form factors formulae the only difference being due to the 
phases P, P'. This difference is quite understandable: in quasi-classical construction we get automatically 
the states which are symmetric with respect to permutation of particles, while in usual form factor formulae 
the states are used which produce S-matrices under these permutations. This difference in normalization 
is responsible for presence of the phases. Notice that the presence of the functions £ in the quasi-classical 
result is due to the contribution of "vacuum" particles. This contribution could not be reproduced in more 
naive approach of the paper || . 

It is clear that we are not very far from the exact quantum formula for the form factors in finite volume. 
The main feature of both quasi-classical formula (|3^ ) and the exact formula in infinite volume is that 
they are given by products of certain multiplier which is independent on particular local operator and finite- 
dimensional integral depending on the local operator. Will this structure will hold for the exact formula in 
finite volume? This is not clear, but this is the only chance for the formula to be efficient. 



8 Appendix. 

In this Appendix we give details of calculation of the quasi-classical limit of Q. Recall that 

ew = e( o>n(i-(A) J ) 

where the zeros ct 2 are real positive. Following Jl7|] introduce 

a(A) = o(A^y 

Using the fundamental equation 

t{\)Q{\) = Q(\q) + QiXq- 1 ) 

one finds: 

loga(A) - lo ga (A< Z - 1 ) - ~ J log(o(/i) + l)dlog Qj^'O (41) 

c 

where the contour C encloses the points aj and <? 2 c| in the plane of /i 2 . 

We know that in the quasi-classical limit the zeros <Tj are dense inside n + 1 intervals corresponding to 
cuts on the fig. 3. Let us present quasi-classically 

o(A) = Oo(A) (1 + Tnrx(X) + 0{Ti 2 )) 
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In the order fi 1 the equations ( fhT| ) give 

^ logao(A) =-^ij log(aofc) + 1) 7^ j ^ 2)2 V (42) 

where the contour C encloses cuts on fig. 3. We know in advance from classical consideration that ao(A) = 
A 2 (A). Let us check, first, that this oq{X) indeed satisfies the equation (|42|). We have 

log(A 2 ( M ) + 1) = ~logA(/x) +logT( M ) 

where T = A + A -1 is entire function without zeros on the cuts. So, the integrals in (E^) with logT(/x) 
disappear leaving an obvious identity. It is clear, on the other hand, that this is the only way to satisfy this 
equation. Notice that the equation does not have solution for arbitrary positions of the branch points A 2 : 
the equations (0) have to be satisfied. 

Consider now the equation ( f4l| ) in the order h 2 . After some simple transformations it can be rewritten 

as 

l-A 2 (Ah 1 d, 1 /"1-A 2 (u) A 

- ■ A- 1 x-(A)^»— logA(A) + - / - ■ ^ gfc)— — Jl (43) 



1 + A 2 (A) v ' dX Tri J l+A 2 (fi) /x 2 - A 2 p 

c+ 

where and the contour C+ encloses all zeros on A 2 + 1 (or T) in the fi 2 plane. This equation shows that 
x(X) allows analytical continuation to the surface E which is realized, as usual, as the plane with cuts (fig. 4). 
Consider the differential 

m !- A2 ( A ) m dA 
^ (A) = ITa^o^T 

The equation (E3) shows that 

1. </?(A) is holomorphic differential on the surface with infinitely many simple poles at the points ±Tj (zeros 
ofA 2 + l). 

2. Since the function x(A) is regular on the first sheet of E (in the upper half-plane) tp(X) has simple zeros 
at the points ifij (zeros of A 2 — 1). 

3. The equation holds: 

ip(X) + tp(—X) = 2id log A(A) (44) 



4. The a-periods of (p vanish. 

Let us show that these four conditions are sufficient to satisfy the equation (^). Consider the canonical 
symmetric second kind differential w(A,/i) = G(A, fi)dXd/i: 

dXd^i 



w(X,fi) =w(p,,fi), w(\,fi)~ — — A~//, / G(A,^)dA = Vj 



It is clear that on our hyper- elliptic surface 

G(A,/0+G(A,-/i)= 2 A/X (45) 
Let us show that the following equation holds for 93 satisfying the conditions 1-4: 

¥3 (A)=i^logA(A) + ^ y <p(ii)Ju>{\,n) (46) 

where C+ (C_ ) enclose zeros on A 2 + 1 in upper (lower) half-planes. Indeed, the contour C+ + C_ encloses 
all the singularities of the integrand except for that at fi = A and that at fi = 00. Obviously residue of the 
integrand at /i = A equals 95(A). Thus using the Riemamm bilinear relation one has 

J <p(n) J w(X, fi) = 2iri<p(X) +lirn R _ >00 J <p(n) f fi) + 

C++C- S + +S- 

+ t ( I w(A, M ) / <Km) - / / w(A, M ) ] 
19 



The latter sum vanishes because the differential ui and ip have vanishing a-periods. The contours S+ (S 1 -) 
are half-circles in the left (right) half-planes of radius R. Using the condition 3 one has 



lim 

R^oa 



cj(A, n) = lim 2 

R — >oo 



S++S. 



2 J idlogA(fi) J 



uj(\,ri = 2nL lim ^^S. 

n^oo dp, 



Recall mow that w(A,/i) is normalized second kind differential with double pole at A = \i while dlogA(A) is 
a normalized second kind differential with singularity at infinity where it behaves as LdX. Hence 

l im ^^M) =L -i dlogA(A) 
m^oo d/i 

which proves the equation (p7j| ). Using (Q) and ( f45| ) and the fact that dlog A(A) does not have singularities 
inside C+ one shows that ( |43| ) follows from (^). Thus tp satisfying the conditions 1-4 satisfies the equation 

©■ ^ 

Let us return to the function x(A). From (M4) one gets 



x(A) + x(-X) 



-2iA^log(A(A)-A- 1 (A)) 



(47) 



Recall that 



Let us present x(A) as 



A(A)-A- 1 (A) = AVP(A)n 1 



3=1 



J*-J. 



x(A) = 2iA 



d / 1 



dX V 2 



log(A/P(A))+logF s (A) 



The function Fs(A) must be regular in the upper half plane, so the simple poles of the RHS of (|47| ) must 
be shared between x(A) and x(— A) in such a way that dlogFs(A) is a differential with simple poles as the 
points —i/A—j with residues 1. Such a differential is not uniquely defined: one can add to it holomorphic 
differentials. These are zero-modes of the equation (fl3|). The origin of these zero-modes is clear: the classical 
surface is parametrized by n moduli Ty, ■ ■ ■ ,T n , the zero- modes correspond to variation of these parameters. 
Since after all we have to impose the Bohr-Sommerfeld quantization condithions (^0|) one can show that the 
final result for Q does not depend on these zero-modes. We require that the ^-periods of dlogFs(A) vanish 
because with this choice the quantum rapidities of solitons don't differ from the classical ones in L — ► oo 
limit. Thus one presents (ilogFs(A) in the form given by the formula (p5|). 
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